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Conventional analysis of in vitro assays of motor proteins rests on the assumption that all proteins 
with the same chemical composition function identically; however molecule-to-molecule variation is 
often seen even in well-controlled experiments. In an effort to obtain a statistically meaningful 
set of time traces that simultaneously avoid any experimental artifacts, we performed quantum- 
dot labeled kinesin experiments on both surface and levitated microtubules. Similar to glassy 
systems, we found that mean velocities of individual kinesin motors vary widely from one motor 
to another, the variation of which is greater than that expected from the stochastic variation of 
stepping times. In the presence of heterogeneity, an ensemble-averaged quantity such as diffusion 
constant or randomness parameter is ill-defined. We propose to analyze heterogeneous data from 
single molecule measurements by decomposing them into homogeneous subensembles. 


Single molecule experiments in the last decades have 
greatly increased our understanding of biomolecules 
through real-time visualization of molecular movement 
m, making it possible to perform time series analy¬ 
sis for individual molecules, as well as to calculate the 
distribution of dynamic variables that had previously 
been difficult to determine in ensemble-averaged mea¬ 
surements. One of the most striking, yet not well ap¬ 
preciated, observations from single molecule measure¬ 
ments is that there are persistent heterogeneities at 
the molecular level in a number of biological systems 
El SHa. In a system with molecular heterogeneity, 
even if all the molecules are chemically identical and 
under the same experimental conditions, a dynamic 
pattern observed in one molecule differs substantially 
from other molecules and this variation is greater than 
the stochastic variation expected from the mean. 

According to the general principle implied by the 
thermodynamic hypothesis in molecular biology (An- 
finsen’s dogma) [8], (i) native states of at least small 
proteins, which can also be extended to RNA El? 
uniquely determined by a given sequence and exter¬ 
nal condition, and (ii) dynamics of biomolecules oc¬ 
curs reversibly with high kinetic accessibility. Both 
conditions are regarded to be critical for biomolecules 
to achieve functional fidelity. In the light of such prin¬ 
ciples, it is rather surprising to find that variation of 
dynamical property is greater than the variation al¬ 
lowed by the stochasticity of rate process, and such 
heterogeneity persists far longer than a biologically vi¬ 
able time scale 013 . Recent single molecule studies, 
such as the docking-undocking transition of ribozymes 
El El and isomerization dynamics of Holliday junctions 


[5], clearly show that dynamics of molecules in action 
can differ drastically from one another. This puzzle 
has recently been illuminated by the variable mean 
velocity of the DNA helicase, RecBCD [6]. Formally, 
kinetics of a single enzyme is stochastic, so it is a priori 
not obvious whether the variation simply reflects such 
stochasticity. Critically, the theoretical analysis to un¬ 
derstand such heterogeneity is poorly developed. Here, 
we discuss when the observed single-molecule variation 
is too large to be generated from molecules function¬ 
ing with the same underlying kinetic parameters. We 
then apply it to kinesin, and discover that such het¬ 
erogeneity exists in the single-molecule function of the 
kinesin motors. 

To preclude the possibility that the observed hetero¬ 
geneity is an undesired outcome of experimental arti¬ 
facts, we performed a few distinct control experiments, 
(i) In bead assays, nonspecific interaction between 
bead and the kinesin tail domain could lead to the 
heterogeneous dynamics. To avoid this possibility we 
attached a bead to the kinesin tail domain via a specific 
streptavidin-biotin interaction (see SI). The trajecto¬ 
ries generated from the bead assays (A^=32, Figj^) 
vividly show velocity variations consistent with motor 
heterogeneity. To determine if the heterogeneous en¬ 
semble of motors could be partitioned into functional 
sub-classes, we developed a new QD/TIRF experimen¬ 
tal setup (Fig.SI), allowing improved (A^=397) statis¬ 
tics. Importantly, in our TIRF-visualized quantum- 
dots (QDs), to which the motor’s tail domain is specif¬ 
ically attached via an antibody, the motors (K560) 
lack the auto-inhibitory portion of the tail [TO], so 
any observed variability does not reflect variable auto- 
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inhibition. Again, velocity variations were observed, 
consistent with motor heterogeneity that is intrinsic to 
the individual molecules, and suggests a need to recon¬ 
sider many facets of the conventional single molecule 
data analysis. 

Kinesin-1 is a motor protein that “walks” proces- 
sively along microtubules by converting chemical po¬ 
tential of ATP hydrolysis into mechanical stepping 
[T]. At saturating ATP condition, it is often stated 
that kinesins move at mean velocity oi V ~ 0.8 
limjsec {=S nm/10 ms) [1]. However, such a descrip¬ 
tion does not make explicit whether the average ve¬ 
locity of 0.8 [imjsec is a common property of all the 
individual kinesin molecules or it is simply a property 
calculated over the heterogeneous population. 

The actual time traces of kinesins (K560) from bead 
assays display substantial heterogeneities in the mean 
velocity (Fig{^). This large dispersion seems unlikely 
to be explained by stochasticity of stepping time com¬ 
prised of binding of ATP, hydrolysis, and ADP re¬ 
lease. While kinesin motors remain on microtubule 
tracks, slow motors are persistently slow and fast mo¬ 
tors are persistently fast. Slow-to-fast or fast-to-slow 
interconversion of velocity is not observed on the time 
scale that a kinesin remains on its track. Since each 
realization of a kinesin time trace is ideally a conse¬ 
quence of the cumulative sum of stochastic yet uncor¬ 
related steps, one can assume that the stepping time 
is drawn from an “independent and identically dis¬ 
tributed” (i.i.d.) random variable. In this case, due 
to the elementary renewal theorem [12], any kinesin 
molecule {a = 1,2, •••A') should have the identical 
mean velocity at long time (or sufficiently large num¬ 
ber of steps), t ^ oc, i.e., limt^oo Vait) = V for any 
a. Evidently, even when the number of steps is large 
enough (t/r ^ O(IO^) > 1), the time traces in Fig{^ 
show persistent and distinct motor movement leading 
to distinct mean velocities, thus violating LLN as in 
glassy systems [iSlIii]. The assumption of i.i.d. is not 
valid among the kinesin time traces. 

To makes this point more straightforward, let us 
consider an ensemble of Poisson walkers that take steps 
with a stepping time distribution '^(t) = 

The positions of walkers along the track are expected 
to evolve as dtx{t) = E + r]{t) where r]{t) is Gaus¬ 
sian white noise obeying P[r]{t)] oc 
with {r]{t)r]{t')) = 2D6{t — P); then the probabil¬ 
ity of finding a walker at position x at time t is 
P[x{t)] = exp [—(x(t) — , and 

hence the mean velocity x{t)/t = V{t) up to time t 
should obey the gaussian-like velocity distribution 


PG[V{t)] 


f * 

{v{t)-vri 
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The variance of V{t) decreases as and the distri¬ 
bution finally converges to S{V (t) —E), suggesting that 
after a sufficient number of steps, any walker ought to 
have an identical mean velocity V. Indeed, an ensem¬ 
ble of simulated time traces with r = 10 ms and d = S 
nm confirms that the mean velocities from an ensem¬ 
ble of Poisson walkers at time t, are distributed as a 
gaussian, centered around V = 0.8 jam/s (Fig.S2a). 
The same conclusion is reached even for a more gen¬ 
eral stepping time distribution, though a single rate 
limiting step ^ e~'^P) is a relatively “worst- 

case” scenario, generating relatively large amounts of 
heterogeneity in velocity. As an example, a step¬ 
ping time r composed of multiple internal substeps 
^ P^e~^P) leads to less dispersed time traces 
(Fig.S2b). As long as the stepping time is i.i.d. and 
has a finite variance {{{St)‘^) < oc), the above conclu¬ 
sion, i.e., P{V{t)) = S[V{t) — F] as t ^ oo is always 
valid, which simply restates the central limit theorem. 

In practice, kinesin motors have a finite proces- 
sivity satisfying exponential travel time distribution, 
Pj^{t) = (F/I/)e“^^/^, as is the case in our simulation 
result where a detachment probability (p = 0.01) is 
imposed on each step (see Figj^). For an ensemble 
of homogeneous motors with exponential travel dis¬ 
tance, the distribution of mean velocities ought to be 
described by incorporating the travel time distribution 
as a weighting factor, which leads to: 


homo 


poo 

{V;V,D,L)= / dtPG[V{t)]Pj;{t) 
Jo 

_ V/L \(V-Vf V 
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^homo(fo is symmetric with respect to V with 

a single peak and power-law tails at both ends. It is of 
particular note that from the scatter plot of (Fq,,I/q,) 
(Figj^) traces with large run length is found predomi¬ 
nantly near the mean velocity, and the traces showing 
the large deviation from V always have a short run 
length. The full width at half maximum (FWHM) 
of Eq§is 2(22/3 - 1)1/2 (v(d/I) 1/2 « 3.1(VD/I)1/2. 
The predicted velocity distribution is indeed realized 
in velocities calculated from simulated Poisson walkers 
(see Fig{^); however, we observe that it fails to ad¬ 
equately describe our experimental data in a number 
of ways. 

Although the violation of LLN is obvious among the 
time traces from the bead assays (Figj^), with only 
A = 32 traces we could not investigate whether there 
might be distinct subgroups of motors with similar 
function. Thus we also used QDs (Fig.Sl) to moni¬ 
tor the movement of individual kinesin motors; this 
experimental geometry not only gave more traces, but 
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FIG. 1: (a) Time traces of single kinesin motors from bead assays (N = 32). The tail domain of individual motor is 

specifically attached to the bead through strepavidin-biotin with 30 % of kinesin/bead binding fraction Scatter plot 
of (Va,Ta) is shown on the right with the corresponding histograms, (b) Computer generated time traces {N = 5000) with 
(r = 10 ms) and detachment probability of p = 0.01 at each step. The histogram of travel distance (right 
panel) is fitted to P{L) = L with L = 0.80 jim (red line); and the velocity distribution (top panel) is to EqMwith 

D = 0.0033 V = 0.81 jj^mls, and L = 0.89 pm (red line), (c) QD-labeled kinesin data on surface-immobilized 

MTs (N = 397). The velocity distribution is fitted to Eq|^ (red line). The green line depicts the cumulative sum of the 
histogram. Meanwhile, the histogram of travel distance (see Eig.S3 for unnormalized one) is fitted from the second bin 
because the trajectories with short run length, which contributes to the first bin, are excluded from the velocity analysis, 
satisfying P{L) = L with L = 0.78 pm. (d) QD-labeled kinesin data on levitated MTs {N = 52). 


was perhaps more natural than the bead/optical trap 
assays that bring motors into contact with the MT 
because individual QD-labeled motors land sponta¬ 
neously on MTs and start to move. 

In qualitative agreement with the bead assay data, 
the time traces and P{V) with multiple peaks from 
TIRE-visualized QD data show the motor heterogene¬ 
ity more convincingly (Fig{^). The time traces in Fig. 

clearly contradict the usual assumption of i.i.d., sug¬ 
gesting that not all the kinesin motors are functioning 
identically, and that the energetic cost of overcoming 
such heterogeneity is not minor. One might wonder 
whether interaction between either the QD or motor 
nearby surface to which the MTs are attached can ac¬ 
count for the observed heterogeneity. If the surface in¬ 
teraction were significant, the velocity of kinesin would 
vary depending which protofilament the kinesin walks 
along. To preclude this possibility, we performed QD- 
assays on MTs levitated above the surface. Evident 
from Fig. we still find the motor heterogeneity 


even in the QD-labeled kinesin time traces from the 
off-surface measurements. The mean velocities of mo¬ 
tors on the levitated MTs (Fig. df) are widespread, 
demonstrating qualitative agreement with those on the 
surface-immobilized MTs. 

Eqj^ that assumes homogeneity of motors, fails to 
explain the overall asymmetric and multiply peaked 
velocity distribution of our QD-labeled kinesin data 
in Fig. [^. Thus, we propose a heterogeneous veloc¬ 
ity distribution by combining homogeneous subpopu¬ 
lations: 

n 

-fhetero(f^) ~ omo {V;Vi,Di,L) (3) 

i=l 

where = 1- The velocity distribution in Figj^ 

can be nicely fit with Eq|^and decomposed into three 
subpopulations (n=3) with (pi = 0.74, p 2 = 0.10, 
03 = 0.16, Vi = 0.85 [imls, V 2 = 0.51 [imjs, 
F 3 = 0.^ pm/ 5 , Di = 0.00^ pmVs, D 2 = 0.0025 
[xw? js, Ds = 0.0043 ixrv? js, L = 0.75 jxm. This sug- 
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gests that there are kinesin subpolulations with fast, 
intermediate, slow mean velocities. The peak around 
is the most dominant, representing ^ 74 % of the 
population. As an alternative analysis method used in 
qunatifying the anomalous diffusion from single par¬ 
ticle tracking na ig, the displacement distribution 
calculated over the ensemble of time traces can also 
reveal the three subpopulations in the ensemble (see 
SI and Fig. S4). 

Although the assumption that all motors move with 
the same mean velocity V at t ^ oo is implicitly made 
in many studies [niisiTn], variability among the mo¬ 
tors obviously exists in the raw data. In the presence 
of motor heterogeneity, the ensemble averaged quanti¬ 
ties, especially diffusion constant {D) and randomness 
parameter (r), are no longer well-defined quantities. 
Here the randomness r is a measure of temporal regu¬ 
larity of motor step [iTl [20] (see SI): 


(t)2 


[xHt)) - {x{t)y 
d{x{t)) 


(4) 


where (...) denotes an average over the ensemble and 
d is the average step size (for kinesin, d ^ S nm). From 
the second expression of Eqj^ it is usually argued that 
{x^{t)) — {x{t))‘^ ^ 2Dt and d{x{t)) ^ dVt as t ^ oc. 
Hence, as t ^ oc, the randomness, r = 2DldV^ is 
ideally a time-independent, dimensionless measure of 
the dynamical fluctuations of the motor, or the disper¬ 
sion of the motors {x{t)) along the track [T7[ [20l |2T]. 
The variation of r, which can be linked to the number 
of rate limiting steps, has been examined as a func¬ 
tion of load and ATP concentration m ns®, and 
it was concluded that kinesin data are best explained 
using a kinetic model with 4 states [22]. However, for 
these analyses based on the randomness parameter to 
be valid, the condition of i.i.d. should be obeyed. 

To further clarify, consider the position of a motor 
protein a that evolves through the following stochastic 
equation dtXa{t) = Ha + where Va is the mean 

velocity, {r]a{t)) = 0 and {r]a{t)r]a{t')) = 2DaS{t - 
t'). The mean square displacement (MSB) at time t 
calculated by averaging over the motor population is 
(see SI): 


{x^it))-{x{t))^ = {{SVf)t^ + 2Dt. (5) 

where {x{t)) = Ea=i ^a(i) = Ea=i ^4* = 
Vt, {{SV)'^) = - F)2 and D is the dif- 

fusion constant averaged over distinct ensembles. Un¬ 
less ((^U)^) vanishes, the quadratic contribution of 
time becomes predominant at t > t* = 2D/{{SV)^). 
Thus, in the presence of motor heterogeneity, the non¬ 
vanishing term renders the randomness pa¬ 

rameter, r, no longer time-independent. Note that the 


Ref. [18] demonstrates the linearity of variance of po¬ 
sition only in the short time limit. 

For an ergodic system, one can calculate the MSB by 
taking either an ensemble average over the time traces 
or a moving time average over a sufficiently long single 
time trace. However, when a system of interest has a 
subsample-to-subsample variation like in glassy mate¬ 
rials, and the system is persistently heterogeneous on 
the physically meaningful time scale (not being able 
to interconvert between different components), the dy¬ 
namics of the system is effectively nonergodic [23] |24] . 
A systematic procedure to analyze data with (weakly) 
broken ergodicity is to identify components (or do¬ 
mains) of the system in its configurational space as 
proposed by R. G. Palmer [5ll23]. 



FIG. 2: (a) Time averaged MSB of QB time traces using 
Eq|^ MSB ~ t for t < 2DalV‘^ and MSB ~ for t > 
2DalV^. (b) Bistribution of = 2DalVa- (c) Scatter 
plot of {Vex,Da) and the corresponding histogram P(U), 
P{D). The data, partitioned into 4 clusters, are depicted 
in different colors. See Fig.S5 for the analysis of the bead 
assay data. 


Eor a motor system with heterogeneity, whose en¬ 
semble properties cannot be represented with a single 
U or B, one can calculate the time averaged MSB 
[24[[25] using 

{{SXa{t))'^)T = j + t) - Xa{t')f dt' 

= + 2Dat ( 6 ) 

with T t for each molecule a, and quantify both the 
directed (V^) and the diffusive (Da) components. The 
log-log plot shows MSB^ t for t < t* and MSB^ for 
t > (Eigj^); a crossover from diffusive to directed 
dynamics is observed at t* = 2Ba/U^ (Eigj^). In ad¬ 
dition to the P{V) in Eig{^, the scatter plot {Va,Da) 
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(Figj^) clearly visualizes the presence of heterogene¬ 
ity with multiple subpopulations. Using a multivariate 
Gaussian mixture model [26] we partitioned {Va^Da) 
into 4 clusters. When the outlier subpopulation (clus¬ 
ter 3) is excluded, the remaining three subpopulations 
show a weak correlation between and The scat¬ 
ter plot of (W: Da) suggests that kinesins with higher 
mean velocity tend to have more variation in stepping, 
reminiscent of the recent finding that the heat released 
during catalytic turnover enhances the enzyme diffu¬ 
sion [27] . 

To summarize, our study provides simple criteria 
to judge the homogeneity of motor time traces: (i) 
Time traces with large run length should have a con¬ 
verging velocity. P{V) ought to be (ii) uni-modal, 
(hi) symmetric with respect to its mean, and (iv) de¬ 
scribed by Phomo(^) (Eqj^ (or should satisfy FWHM 
^ 3 . 1 {VD/L)^Dy If a distribution satisfies the cri¬ 
teria (ii) and (iii), yet heterogeneity is still suspected 
then one can check the criterion (iv) (Statistical assess¬ 
ment of a possible difference between the two distribu¬ 
tions, P{V) and Phomo(E) can be made, for example, 
using the Kolmogorov-Simirnov test). Violation of the 
above criteria most likely indicates the presence of mo¬ 
tor heterogeneity. Once heterogeneity is identified, we 
propose to decompose the heterogeneous ensemble of 
kinesin motors into a finite number of homogeneous 
sub-ensembles. 

Care should be taken in developing theories for het¬ 
erogenous samples, as has been pointed out in the 
problem of phase separations of practical polymer 
samples with polydispersity [28l [29]. In biology, al¬ 
though importance of heterogeneity is increasingly ap¬ 
preciated at the cellular level [30l [31] , it is surprising 
to observe such characteristics at the level of single 
biomolecules. The structural origin of functional het¬ 
erogeneity or multiple native states mi?] is not clear. 
Yet, all the data from different control experiments dis¬ 
cussed here unambiguously indicate that kinesins can 
have multiple functional states with distinct capability 
of processing ATP, which gives rise to fast, intermedi¬ 
ate, and slow groups of motors. It is plausible that 
while the conformational cycle of kinesin is “driven” 
by chemical potential from ATP hydrolysis, functional 
states of motile kinesins are dynamically pinned or 
separated by large free energy barrier over which no 
thermodynamic path can easily connect one state with 
another. Although it is not easy to test this hypoth¬ 
esis for kinesins due to the short processivity, recent 
studies of Holliday junctions [5] and RecBCD [6] have 
shown that subensemble-to-subensemble interconver¬ 
sion can be induced by depleting cofactors (Mg^+ and 
ATP, respectively) for a finite amount of time. 

Regardless of the cause of the heterogeneity, it is im¬ 


portant to appropriately characterise it, because one 
goal of single molecule studies is to measure proper¬ 
ties to allow prediction/calculation of ensemble func¬ 
tion. Ensembles of heterogeneous molecules are likely 
to function quite differently from ensembles of homo¬ 
geneous ones. For instance, heterogeneous motors may 
interfere more with each other, and in general systems 
combining multiple motors with different velocities, 
unexpected ensemble properties can emerge [321135] . 
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SUPPLEMENTARY INFORMATION 

Experimental Methods. 

Protein: To ensure that the population of kinesins 
are chemically “identical”, we expressed the func¬ 
tional, truncated kinesin (K560) [36] in terrific broth 
and purified it as described in Ref. [36l [37] . Note 
that the lack of light chains (tails) prevents potential 
interactions between light and heavy chains that may 
alter function. Further, we minimized the chance of 
phosphorylation, another potential source of chemical 
heterogenity, by expressing the kinesins in E. coli. 
Furthermore, to select only the functional kinesins, 
the stock protein was selectively purified via MT 
binding and release in the presence of AMP-PNP (the 
nonhydrolizable ATP analogue). 

Quantum Dot (QD) Motility Experiment: Trun¬ 
cated kinesin-1 was specifically recruited to quantum 
dots via its genetically encoded C-terminal His-tag 
[38] . To achieve the specific linking of kinesin to 
the cargo, streptavidin quantum dots (QD-655- 
streptavidin conjugate. Life Technologies) were 
labeled with biotin conjugated pent a-his antibody 
(Qiagen) in molar ratio of 1:1::AB:QD. The incu¬ 
bation temperature of antibody with QDs was 4C 
and lasted for 1 hour. The AB labeled QD surface 
was blocked with 4mg/ml casein (sigma Aldrich, 
C8654-500G) in the motility buffer (80 mM Pipes pH 
6.9, 50 mM CH3C02K, 4 mM MgS 04 , 1 mM DTT, 
1 mM EGTA, 10 jjM taxol, 4mg/ml casein) for Ih at 
4C to reduce nonspecific binding. With this blocking 
procedure K560 binding to casein coated streptavidin 
QDs alone (without antibody) was negligible (reduced 
by 50 times) compared to the same for AB coated 
QDs. K560 with His tag on its truncated tail was 
incubated with anti-His antibody tagged QDs in the 
molar ratio of 1:40 at RT for 12 minutes. Before using 
the mixture to test motility it was supplemented 
with 1 mM ATP and oxygen-scavenging system. 
The high ratio of K560:QD (=1:40) was chosen to 
ensure the likelihood of recruiting a single kinesin 
to the QDs [37]. A sample chamber assembled with 
taxol stabilized microtubules made from bovine 
brain tubulin (cytoskeleton) was used for motility 
experiments. Consistent with the hypothesis that 
individual QDs were at most attached to a single 
kinesin, the distribution of run-lengths was well fit 
by a single decaying exponential distribution, with a 
mean travel of 0.78 microns (Fig. S2). This travel 
distance is what would be expected for single-motor 
travel. 

TIRE Imaging and Traeking: The QDs were imaged 


via custom built objective based total internal reflec¬ 
tion fluorescence microscope (Nikon 1.49NA, lOOX, 
and 488 nm laser Ti:Sapphire, Coherent). The time 
lapse movies were recorded at 31fps using EMCCD 
camera (Phot metrics Quant EM 512SC). Tracking 
analysis was carried out using a custom matlab 
program (Gross Lab) that identified QD positions in 
the images via 2D Gaussian fitting of their intensity 
profile. 

Levitated mierotubule assays: To eliminate concerns 
about surface interactions, we developed elevated mi¬ 
crotubule assay. In this assay 500nm carboxyl ter¬ 
minated polystyrene beads were incubated with enzy¬ 
matically dead kinesin (a mutant, E237A in hKIF5A) 
at RT for 5 minutes. The single E237A mutation in 
the enzymes switch 2 region precludes ATP hydrolysis, 
rendering the motor immobile even in the presence of 
saturating ATP concentrations thus acting like a MT 
anchor. After recruitment the beads surface was then 
blocked with casein, and the beads were spun down at 
5000g and re-suspended in casein buffer, to eliminate 
free kinesin in from buffer. These beads were then 
flown into a flow cell, and allowed to stick to the pre¬ 
cleaned poly lysine coated coverslips for 6 minutes. At 
the end of 6 minutes the new buffer with high casein 
levels (5.5 mg/ml) was flowed in, to reduce nonspecific 
binding of MTs and QDs to coverslip surface. Gasein 
coating made the surface minimally reactive, so that 
when taxol-stabilized microtubules were flown into the 
coverslip, they would not stick to the coverslip, and 
could easily be washed out if the beads with mutant 
kinesin were not present. However, in the presence 
of mutant kinesin coated beads, taxol-stabilized MTs 
ended up being captured by the mutant kinesins on 
the bead, so that the MTs were suspended between 
beads, but above the surface. As this was a stochastic 
process, some MTs ended up close to the surface, and 
some further away. 

Once the microtubules were in place, we now 
flowed in QD-labeled kinesin motors (as above) 
and recorded their motion via wide-held semi-TIRF 
imaging, as shown in Fig. Hit- Importantly, since 
the MTs were at different heights above the surface, 
to measure the height (distance from surface) of a 
particular QD, we used the elipticity in the point 
spread function of the QDs (FigjS^). In brief, by 
putting in a cylindrical lens (f= 1000mm, at 80mm 
from GGD chip) in front of the camera, we could 
induce astigmatism similar to 3 dimensional STORM 
imaging [39]; by quantifying the amount of distortion 
for a given QD, we could determine its height. To 
calibrate this system, the detector was focused on the 
coverslip surface, and then the surface was moved in 


known increments in the Z direction (via a piezo- 
controlled stage), and we quantified how the images 
of surface-fixed QDs changed with Z displacement 
(Fig. [ST]^). For our studies the elipticity of the “off- 
surface” QDs reflects an average QD height of 250 nm. 

Analysis of kinesin time traces using displace¬ 
ment distribution. The displacement probability 
distribution, a quantity often calculated for the 
purpose of revealing the anomalous diffusion in single 
particle tracking, could be used as an alternative 
method for analyzing our in vitro kinesin data. The 
relative increment of displacement measured with a 
uniform time interval of 60 milisecond calculated for 
the ensemble of our QD-labeled kinesins on surface- 
immobilized MTs shows non-gaussian distribution, 
fitted to an exponential for large displacements 
(see Fig]S4[). The exponential tail is similar to the 
data analyzed for the colloidal beads diffusing on 
lipid tubes or particles diffusing in F-actin network 
in Ref. [15]. What is even more interesting is the 
displacement distributions calculated at larger time 
intervals (t = 300 ms and 600 ms. See Fig|S4[), which 
apparently capture the presence of the three distinct 
components in the ensemble of the kinesin time 
traces. Since the displacement at a large time interval 
effectively becomes equivalent to the global velocity 
measured for individual motors, the three peaks 


observed in G{Sx, t) at large t indicate that there 
are three distinct subpopulations in the ensemble of 
kinesin motors. 

Meaning of the randomness paramter, r. 

From the definition of the randomness parameter r 


and {r'^) = an exponential poisson 

process '^(t) = leads to r = 1; and 

'0(t) = ln\ to r = l/(n + 1). If r ^ 0, 

there is little variation in each dwell time; thus the 
motor step is clock-like. By contrast, if r ^ 1, the 
duration of individual steps is highly irregular. 

Derivation of Eq.5. Consider a motor protein a 
that evolves through the stochastic equation dtXa{t) = 
Va + Va{t)^ where W is the mean velocity, r]a{t) is 
the Gaussin random noise with zero mean. From 
the formal solution of Xa{t) = Vat + /q drrjai^) with 
^a(O) = 0? fhe ensemble averaged displacement is 
{x{t)) = Vt where {x{t)) = is the en¬ 

semble average displacement of motors at time t, and 
V = ^ ^a=i Then, the mean square displace¬ 
ment (MSD) for the motor ensemble is: 


{x^t)) - {x{t)f 




^^iVa-vrr + 

a=l 

{{5Vf)r + 2Dt. 


1 y ^ r r 

n'Z dT2ria{Tl)Va{T2) 

^ a=l "'0 -'0 


(S2) 


where we define the variance of the mean velocities as {{6V)‘^) = Y^a=i(^<^ ~ Tor the noise related term 

from the second to the third line of Eq|S2[ we use the property of Gaussian random noise, whose correlation 
satisfies ^ = ‘2Dp6{ri — r 2 ), where the subscript p denotes the index of a subensemble {p = 

1 ,..., He); hence. 


nt nt nt ^ 

VE / / dT2ria{Ti)ria{T2) = —Yi / (Ti)r?i (tz) 

^ a=l do Jo ne Jo Jo 

2 pt nt 

= — Yi / dx22Dp5{Ti - T 2 ) 

Jo Jo 

= 2Dt (S3) 

where X]p=i 1 = and D = y X]p=i is the diffusion constant averaged over the distinct subensembles. 
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FIG. SI: Geometry and calibration of off-surface measurements, (a) Diagram of experimental apparatus, showing 
location of cylindrical lens (GL) in front of camera, to induce distortion when the QD is out of focus, (b) Diagram of 
experimental geometry. Half-micron beads coated with mutant kinesin (which rigor-binds to MTs) are attached to the 
coverslip, and microtubules are subsequently flowed in, and stick to the beads, ending up suspended between beads above 
the surface, (c) Quantification of asymmetry in QD image, as a function of the QDs distance from the plane of focus. 
The extent of asymmetry was used to detect QDs moving on MTs either close to or far-from the surface. The error bars 
are SEM, estimated by tracking the position and intensity profiles of 20 QDs in the time lapse images recorded during 
piezo Z motion. 
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FIG. S2: Gomputer-generated time traces (iV = 1000) by using stepping time distributions: (a) — r and (b) 

where ts = t/5 with r = 10 ms and step size of d = 8 nm. The histogram of mean velocity 
(inset of (a)) is fitted to Eq.l in the main text with V = 0.8 iimjs and D — 0.0032 [ivn? js at t = 10 sec. P{V)s at t = 2 
and 50 sec are also plotted. 
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FIG. S3: Histogram (unnormalized) of travel distances 
from QD assays. The decay length from the single¬ 
exponential fit was L = 0.78 /xm. 



6x (nm) 

FIG. S4: Displacement distribution G{Sx,t) calculated 
over the ensemble of kinesin time traces for three differ¬ 
ent time interval t. The tail region {5x ^ 1) of G{6x, t = 
60 ms) is described with an exponential function. 
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FIG. S5: MSD analysis for bead assays with specific attach¬ 
ment. (a) MSD of time traces from temporal moving aver¬ 
age defined in Eq.6. MSD versus time in logarithmic scale. 
MSD ~ t for t < 2DalV^ and MSD ~ for t > 2DcxlV^. 
(b) Distribution of the cross-over time = 211^/1/^. (c) 
Scatter plot of (Va,Da) and the corresponding histogram 
P{V), P{D). 


















